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Abstract 

We investigate the fraction of metal nuclei in the relativistic jets of gamma-ray bursts associated 
with core-collapse supernovae. We simulate the fallback in jet-induced explosions with two-dimensional 
relativistic hydrodynamics calculations and the jet acceleration with steady, radial, relativistic magnetohy¬ 
drodynamics calculations, and derive detail nuclear composition of the jet by postprocessing calculation. 

We found that if the temperature at the jet launch site is above 4.7 x 10®K, quasi-statistical equilibrium 
(QSE) is established and heavy nuclei are dissociated to light particles such as ^He during the acceleration 
of the jets. The criterion for the survival of metal nuclei is written in terms of the isotropic jet luminosity 
as < 3.9 X 10^°(i?i/10'^cm)^(l -\-ai) erg s“^, where R, and a-, are the initial radius of the jets and the 
initial magnetization parameter, respectively. If the jet is initially dominated by radiation field (i.e., a\ <C 1) 
and the isotropic luminosity is relatively high (Tj®° ^ 4 x 10^^ ergs“^), the metal nuclei cannot survive in 
the jet. On the other hand, if the jet is mainly accelerated by magnetic field (i.e., cti ^ 1), metal nuclei 
initially contained in the jet can survive without serious dissociation even for the case of high luminosity 
jet. If the jet contains metal nuclei, the dominant nuclei are ^®Si, and and the mean mass number 
can be (A) ~ 25. 

Key words: gamma-ray burst: general — magnetohydrodynamics (MHD) — nuclear reactions, nu¬ 
cleosynthesis, abundances 


1. INTRODUCTION 

Gamma-ray bursts (GRBs) are one of the most ener¬ 
getic phenomena in the universe. They radiate enormous 
energies of the order of 10®^ ergs mainly in the form of 
gamma-rays with the short duration, typically of ~ 10 s. 
Although the detailed radiation mechanism of gamma- 
rays is under debate, the radiation is thought to be origi¬ 
nated from the ultra-relativistic jets with T ^ 100, where 
r is the Lorentz factor of the jets (e.g., Meszaros 2006). 
Furthermore, it is known that GRBs with duration longer 
than 2 s, called long GRBs, are associated with highly 
energetic type Ic supernovae (e.g., Iwamoto et al. 1998). 
This implies that the relativistic jets are launched at deep 
inside of massive progenitor stars after the onset of the 
core-collapse. 

Nuclear composition of GRB jet have been studied by 
several authors. Pruet et al. (2002), Lemoine (2002) and 
Beloborodov (2003) investigated within the framework of 
the standard fireball model (e.g., Meszaros 2006), assum¬ 
ing that GRB jet initially consists of free nucleons. The 
nucleons can recombine into deuterium and/or a parti¬ 
cles as the jet expands and cools but metal nuclei heavier 
than carbon can not be produced. Thus, they conclude 
that the GRB jet consists only of light nuclei. 

However, there is a caveat that the initial composition 
of the jet is not necessarily dominated by the free nucle¬ 
ons because the mechanism to launch a well-collimated 


relativistic jet is not specified. The proposed mecha¬ 
nisms include neutrino annihilation (e.g., Woosley 1993; 
MacFadyen & Woosley 1999) and magnetic field (e.g., 
Blandford & Znajek 1977; Brown et al. 2000; McKinney 
2006). In fact, some observations of GRBs suggest that 
their relativistic jets were initially dominated by the mag¬ 
netic field energy flux (e.g., Zhang & Pe’er 2009; Guiriec 
et al. 2011). When the relativistic jets are launched by the 
magnetic held, the jets could involve heavy nuclei such as 
®®Ni (Fujimoto et al. 2008) and thus the nuclear compo¬ 
sition of the jets may be different from the one that is 
derived within the framework of standard hreball model. 

Horiuchi et al. (2012) investigated survival of metal nu¬ 
clei in relativistic jets. They analytically estimated the 
conditions to survive photodisintegration and spallation 
at the base of the jets, in the accelerating jets, and at 
the emission region of GRB. They argued that the nuclei 
can avoid the destructions for a range of jet parameters. 
For examples, the metal nuclei can survive at the base 
of the jet if the radius of the central engine is greater 
than ~ 10® cm or if the radiation luminosity is less than 
~ 10"^® erg/s. They also investigated the possibility for en¬ 
trainment of metal nuclei from stellar material during the 
jet propagation in the progenitor star and suggested that 
the entrainment is possible depending on the model pa¬ 
rameters. However, they did not calculate detail nuclear 
reactions and did not present final nuclear compositions. 

In this paper, we calculate non-equilibrium nuclear re- 
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actions during the fallback and jet acceleration, based on 
the thermal histories derived by the relativistic hydro¬ 
dynamics calculations and the steady, radial, relativistic 
magnetohydrodynamics calculation. Then we present the 
criterion for the metal nuclei to survive without serious 
dissociation. For the initial composition, we assume that 
the GRB jet is initially made from falling matter during 
a relativistic jet-induced explosion. We adopt Wolf-Rayet 
stars proposed to be progenitors of GRB-SN and circum¬ 
stantially treat the jet acceleration by thermal and/or 
magnetic pressure gradient and nucleosynthesis in rela¬ 
tivistic jets. 

This paper is organized as follows. In Section 2, we de¬ 
scribe the method and the model of GRB jet. In Section 
3, we present fallback in jet-induced explosions, hydro- 
dynamic properties of the accelerating jet, and the final 
nuclear composition of the jet. The conclusion is pre¬ 
sented in Section 4. Finally, the discussions are presented 
in Section 5. 

2. METHOD MODEL 

2.1. Method 

We calculate the nuclear composition of magnetized 
GRB jets with following two steps: (1) We follow fall¬ 
back in a relativistic jet-induced explosion with a two- 
dimensional relativistic hydrodynamics calculation in or¬ 
der to derive the initial composition of jets. (2) We follow 
the acceleration of magnetized GRB jets after the launch 
of jets with a steady, radial, relativistic magnetohydrody¬ 
namics calculation, assuming an interaction between the 
jet and stellar mantle does not influence on the accelera¬ 
tion of the GRB outflow. 

According to thermodynamical histories during the fall¬ 
back and acceleration, nuclear reactions are calculated 
as a postprocessing (e.g., Hix & Thielemann 1996; Hix 
& Thielemann 1999). The nuclear reaction network in¬ 
cludes 281 isotopes up to ^®Br. Here, we adopt an equa¬ 
tion of state for relativistic particles, p = eint/3, which 
means 7 a =4/3 where 7 a is the adiabatic index, and as¬ 
sume that the internal energy is dominated by contribu¬ 
tion from photons and pairs. Temperature T is derived 
by following equation (e.g., Freiburghaus et al. 1999): 

ei„t = aT4(l + -^^^), (1) 

where a = 7.57 x 10“^® erg cm“^ K~'^ is the radiation 
constant and Tg = T/(10® K). We note that the non- 
relativistic gas pressure is negligible compared with radi¬ 
ation pressure because of the high entropy. 

2.1.1. Fallback in relativistic jet-induced explosions 

In jet-induced explosions, a considerable fallback takes 
place along the equatorial plane (e.g., Maeda & Nomoto 
2003) and thus an interaction between the jets and co¬ 
coon and the stellar mantle determines which mass ele¬ 
ments fallback to a central remnant. Therefore, a numer¬ 
ical simulation is required to correctly treat fallback in 
explosions with relativistic jets. Hence, we calculate rela¬ 
tivistic jet-induced explosions of C-l-0 stars with the use 
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of a two-dimensional relativistic Eulerian hydrodynamic 
code with the Newtonian self-gravity (Tominaga et al. 
2007; Tominaga 2009). 

Since the GRB jets consist of the falling matter, the 
initial composition of the outflow is set to be an integra¬ 
tion of falling matter. However, as the freefall time of the 
matter at the inner boundary is short (< 0.1 s), materials, 
that fall through the inner boundary well before the initi¬ 
ation of the jet injection, are likely to have been accreted 
to the central remnant before the jet injection and not 
to be re-ejected as the relativistic jets. Therefore, in this 
paper, only the matter falling after the initiation of the 
jet injection are assumed to be re-ejected and integrated 
as the initial composition of the outflow. 

As the material falls, the temperature increases and nu¬ 
cleosynthesis may take place. Nucleosynthesis during the 
infall from the presupernova location to r = Ri is calcu¬ 
lated with a thermodynamical history taking into account 
heating due to the infall and cooling due to the Urea pro¬ 
cess (Bisnovatyi-Kogan 2002). 

2.1.2. Steady relativistic magnetized outflow 

We treat a GRB jet as a steady, radial, magnetized 
outflow with efficient magnetic dissipation. We assume 
that magnetic fields in the outflow are dominated by a 
toroidal component and that the field efficiently dissipates 
via magnetic reconnection as the outflow expands. Such 
efficient dissipation of magnetic fields creates strong mag¬ 
netic pressure gradient which enables a direct conversion 
of magnetic field energy into kinetic energy (Drenkhahn 
2002; Drenkhahn & Spruit 2002). We neglect gravita¬ 
tional force because the gravitational energy is small com¬ 
pared with the radiation energy or magnetic field energy 
in the models considered here. 

The mass, momentum, and energy conservation equa¬ 
tions for the steady, radial outflow with a toroidal mag¬ 
netic field are described as follows (e.g. Lyutikov & 
Blandford 2003): 

= 0 , ( 2 ) 
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where r, p, F, v, p, w, and b are distance from the center, 
proper mass density, outflow Lorentz factor, outflow ve¬ 
locity, pressure, enthalpy, and toroidal component of the 
magnetic four-vector, respectively. The enthalpy is de¬ 
fined hy w = p-\- Cint + P where e-mt is the internal energy. 
Here we adopt the unit system in which the speed of light 
is unity. 

We employ a following evolution equation for the mag¬ 
netic field which includes dissipation of non-axisymmetric 
magnetic field produced by inclined rotator (Drenkhahn 
2002 ): 
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where r<iis is the dissipation time scale derived as follows: 
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where ua = bj y/w is Alfven four-velocity, 17 is an angular 
frequency of central object, and e is a dimensionless factor. 

From Equation (2)-(6), the evolution of the Lorentz fac¬ 
tor can be written as 

i9r _ 7a-c^r3 j'2p 

dr (r2-r^)(ui- 7 ap) 1 r "faVTdis j’ 

where Ff is a Lorentz factor corresponding to phase ve¬ 
locity of a fast magnetosonic wave and can be expressed 
as (e.g., Lyutikov & Blandford 2003) 
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Equation (7) indicates that the flow is accelerated only 
when F > Ff and the acceleration of the flow is infinity 
when F = Ff. These behaviors of the flow is due to the 
absence of the gravity in the above formulation. Thus, 
although the magnetic field possesses a part of the total 
energy, we assume that the flow is accelerated as if there 
is no magnetic field and the magnetization parameter de¬ 
fined by cr = b'^/w is constant until F reaches 2rf. Then, 
the flow is accelerated with Equation (7) at F < Ff. 


2.3. Model 


We parameterize the GRB outflow with six parameters: 
isotropic energy deposition rate Lf°, initial radius of the 
outflow i?i, initial Lorentz factor F;, maximum Lorentz 
factor Fuiax, angular frequency of central object (includ¬ 
ing dimensionless factor) ef2, and initial magnetization pa¬ 
rameter (Ti = bf /wi, where 6; and Wi is the initial troidal 
component of the magnetic four-vector and the initial en¬ 
thalpy, respectively. 

In this paper, we investigate the dependence of final 
nuclear composition on i?i, and Ci for the C-fO star 
models with metallicity Z = 0 and 0.02. We adopt progen¬ 
itor stars constructed from C-l-0 cores of 40Mq stars with 
Z = 0 and Z = 0.02 (Umeda & Nomoto 2005) by attaching 
C-fO envelopes in hydrostatic and thermal equilibrium. 
The envelope connects to the core structures continuously 
and smoothly in the first order differentials (e.g., Saio et 
al. 1988) and extends down to a density p= 10~^° g cm“^. 

The density at the O layer is one order of magnitude 
lower in the model with Z = 0.02 than in the model with 
Z = 0. The O layer can be divided into two layers accord¬ 
ing to the abundance of C and Mg, 0-|-Mg and O-fC lay¬ 
ers. The boundaries between the two layers for the models 
with Z = 0 and Z = 0.02 are 3 x 10® cm and 8 x 10® cm, 
respectively. 

The parameter ranges of the outflow are as follows; (1) 
We adopt Lf °2 = Lf°/{10^'^ erg s-^) = (4, 10) because the 
GRB isotropic luminosity is normally-distributed in the 
range of L]^° ^ xo5i-54 gj.g g-i (Ghirlanda et al. 2010). 
(2) Ri corresponds to the size of central engine that is 
constrained by time variability of GRB prompt emission. 
Since an interval of the time variability At is of the order 


of 10“® s, Ri is presumed to be smaller than cAt ~ 10® cm. 
Therefore, we vary Ri for a range of Ri^r = Ri/(10^ cm) = 
1 ~ 10. (3) (Tj is determined by the mechanism to launch 
relativistic jets. We construct models with cTj = 0, which 
corresponds to the standard fireball model, and with a 
range of Ui = 0.1 ~ 10 , which corresponds to magnetized 
jet models. We note that the strength of magnetic held 
realizing adopted Ci is ~ 10^®“^® Gauss at Ri. (4) For 
the other parameters, we hx F„iax = 100 referring to a 
requirement to avoid compactness problem (Piran 2004), 
and expediently set ell = 10 ® s“® assuming e = 0.1 and 
fl = 10^ s“®. We note that efl is rather uncertain. Initial 
Lorentz factor of the outflow is hxed to be Fj = 1.23, which 
corresponds to the sound velocity of ultra-relativistic flu¬ 
ids Cs = l/-\/3. We name the model with (^“ 52 , ^)=(4, 
0) Model A, (L!-, Z)=(10, 0) Model B, (Lj- , Z)={A, 
0.02) Model C, and Z)=(10, 0.02) Model D. The 

model parameters are summarized in Table 1. 

In the two-dimensional relativistic hydrodynamics cal¬ 
culation, we set parameters of the relativistic jets corre¬ 
sponding to the parameters of GRB outhow. The rela¬ 
tivistic jets are injected from the inner boundary at an 
enclosed mass 1 . 4 M 0 corresponding to a radius Rin = 
900 km. The initial half angle of jet is set to 15° (Zhang 
et al. 2003) and thus energy deposition rates of jets Lj is 
set to be = Lj/lO®^ erg s“® = 1.5 and 3.4, which give 
Lj ®52 = 4 and 10, respectively. We employ the Lorentz fac¬ 
tor of the jet at the inner boundary as F = 2.5 by reference 
to § 3.2 and the energy density of the jet as to accelerate 
the jet to F = 100 . 

3. RESULTS 

3.1. Initial composition 

Figures la and lb show the falling regions of the models. 
The matter at the intersection between the circumference 
of the falling region and the jet axis is located at the inner 
boundary when the jet injection initiates. In the models 
with Z = 0 (Model A and B), the jet injection initiates 
at early time for the model with = 10 (Model B), 
while it has to wait until the large portion of the material 
in the O layer falls in the model with Lj ’^52 = 4 (Model 
A) because the ram pressure of the jets cannot overcome 
that of the falling matter in the O layer (see Maeda & 
Tominaga 2009). As a result, the falling region of the 
model with Lj '^52 = 4 is more extended than that of the 
model with Lj ’*52 = 10. On the other hand, in the models 
with Z = 0.02 (Model G and D), due to the low-density O 
layer the jet injections in both of the models with Lj ®52 = 4 
and 10 initiate at similar epochs and thus their falling 
regions are similar. The initial composition of the GRB 
outflow is derived from the integration of matter at the 
filled region because the matter at the shaded region is 
likely to be accreted to the central remnant before the 
initiation of the jet injection. 

The composition of a mass element falling to r = Ri 
depends on the maximum temperature that is higher for 
smaller Ri and for the matter locating initially at the in- 








Table 1. Summary of the model parameters 


Model 

L!™ 

(10®^ erg s“^) 

(10®^ erg s“^) 

Z 

r 

^ max 

Ti 

efi 

(s-i) 

A 

4 

1.5 

0 

100 

1.23 

10^ 

B 

10 

3.4 

0 

100 

1.23 

10^ 

C 

4 

1.5 

0.02 

100 

1.23 

10^ 

D 

10 

3.4 

0.02 

100 

1.23 

10^ 


z=o (a) 

Model A 
Model B 

0+C 



0 2 4 6 8 10 



0 2 4 6 8 10 

R/10® [cm] 


Fig. 1. Initial locations of the falling mass elements (filled 
and shaded regions), for (a) Z = 0 models (Model A and B) 
and (b) Z = 0.02 models (Model C and D). The shaded and 
filled regions represent the mass elements fell before and af¬ 
ter the initiation of the jet injection, respectively. We as¬ 
sume that the outflow consists of the materials initially at 
the filled regions. The color of the regions represent the 
models with ^ (green) and — 10 (red) and the 

background circles represent the boundaries between an in¬ 
ner O+Mg layer with JA(^^Mg) > A(^^C) and an outer O+C 
layer with ^(^^Mg) < 



Fig. 2. Mean mass number of initial composition of the out¬ 
flow. 


ner layers. The mean mass numbers (A) of the initial 
composition of the outflow are shown as a function of 
i?i in Figure 2. The Model B has small (A) = 1 — 6 for 
i?i 7 = 1 — 10 because the jets consist of the matter ini¬ 
tially locating in the inner layer, while the Model A can 
have (A) > 15 for i?i _7 > 5.5 because the heavy nuclei can 
survive in the outer matter that falls to the outer Ri. The 
most abundant nucleus in the models with (A) > 15 is ^®Si 
due to the O burning during the fallback. The turnover at 
i?i ^7 ~ 8 stems from the fact that the maximum temper¬ 
ature of the mass elements in the outer layer is not high 
enough to ignite 

On the other hand, the jet injection in the models with 
Z = 0.02 (Model C and D) are initiated at early time and 
the outflow contains the matter initially locating in the 
O+Mg layer. The falling region are similar in both models 
and thus resultant (A) of the models are similar. In the 
models with Z = 0.02, the presupernova temperature of 
the infalling matter is lower than in the models with Z = 0 
due to the low-density O layer. Therefore, (A) > 15 is 
realized at Ri > 7.3 and the most abundant nucleus is ^®Si 
in these models. 

3.2. Outflow dynamics 

Hydrodynamical properties of the outflow for the mod¬ 
els with = (4,6,0) and (4,6,3) are shown in 

Figures 3(aj-3(d). 

Figures 3(a) and (b) show evolution of the Lorentz fac¬ 
tor as functions of r and fluid proper time tp which is cal¬ 
culated by an integration tp = l/(vr) dr, respectively. 
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Radius r [cm] Fluid proper time tp [sj 


Fig. 3. Hydrodynamical properties of the models with 
(L!®g 2 i ^i,7i O'!) = (4,6,0) and (4,6,3): (a) radial evolution of 
Lorentz factor, (b) time evolution of Lorentz factor, (c) radial 
evolution of the temperature, and (d) time evolution of the 
temperature. 

In the model with (Tj = 0, the Lorentz factor evolves lin¬ 
early with radius until it reaches Lmax as expected from 
the standard fireball model. The acceleration takes place 
at tp ^ 10“^ — 10“^ s. On the other hand, in the models 
with (Ti = 3, the outflow is accelerated initially by thermal 
pressure like the fireball model and later by magnetic field 
with the time scale of Tdis ■ The Lorentz factor evolves 
more slowly than the model with Ui = 0 . 

Figures 3 (c) and (d) show evolution of the tempera¬ 
ture as functions of r and tp, respectively. In the both 
models with Ui = 0 and 3, the temperature decreases ex¬ 
ponentially with tp, during the phase in which the Lorentz 
factor evolves linearly with radius. This stems from the 
relations, T oc r, T oc r“^, and dtp = dt/T ~ dr/cT. While 
the temperature in the model with cri = 3 is lower than 
that in the model with crj = 0 before the acceleration, it 
decreases more slowly and becomes higher at r > 10 ^° cm 
and tp > 10“^ s than in the model with o-j = 0. This is 
because the acceleration is slow and the magnetic field en¬ 
ergy is converted not only to the kinetic energy but also 
to the thermal energy in the model with cri = 3. 

3.3. Nuclear composition 

Figures 4 (a) and (b) show time evolution of nu¬ 
clear composition of the outflow for models with 
.Ri, 7 ,CTi) = (4,0,6,0) and (4,0,6,3), respectively. 
Nuclear reaction ceases at tp ~ 10“^ sec in both models 
(a) and (b) because the temperature after the epoch falls 
below ~ 10® K (Figure 3d). The epoch corresponds to 
r~ 10^° cm for the model with (Lj®52,^,.Ri,7,t7i) = (4,0,6,0) 
and r ~ 4 X 10® cm for the model with (Lj™ 2 ,^,.Ri, 7 ,o’i) = 
(4,0,6,3). 

In the model with (Lj® 52 ,-^,.Ri, 7 ,o'i) = (4,0,6,0), which 
corresponds to the standard fireball model, heavy and in¬ 
termediate mass nuclei are almost dissociated until tp ^ 
10 “® s and ^He and d are synthesized at tp > 10 “® s as 
the outflow cools. Only light nuclei and free nucleons, ^He, 
d, p, and n, remain in the outflow after the acceleration 



Fig. 4. Time evolution of nuclear composition for 
the models with (a)(Lj®g 2 ;-^i ,75 ^i) — (4, 0, 6, 0) and 
(b)(Li^g2,Z,/?i,7,cri) = (4,6,6,3). 


(Figure 4a). This demonstrates that the standard fireball 
model with (^.^ 52 ,2', Ui) = (4,0,6,0) destroys heavy 
nuclei even if they are initially contained in the outflow. 

On the other hand, in the model with (Lj®52,^,.Ri,7,<7i) = 
(4,0,6,3), whereas some nuclei such as ®®Ni are dissociated 
at tp < 10 “® s and re-synthesized at tp > 10 “® s, other 
nuclei such as ®®Si survive without dissociation and remain 
abundant after acceleration (Figure 4b). This illustrates 
that the metal nuclei can survive in the magnetized jet. 

Figures 5(a)-5(d) show mean mass number of final nu¬ 
clear composition as functions of i?i and crj for the model 
A, B, C, and D (see Table 1). 

When a condition 


Tg > 3.6 



1 

33.3 


(9) 


is satisfied, quasi-statistical equilibrium (QSE) is attained 
(e.g., Woosley et al. 1973). The white solid line in 
Figures 5(a)-5(d) represents the criterion to establish 
QSE. Therefore, in the models with lower Ui and smaller 
i?i than the white solid line, the composition can be de¬ 
scribed by QSE and the metal nuclei are almost destroyed 
due to the high entropy. This criterion is well fitted by a 
contour of initial temperature Ti^g = Ti/IO^K = 4.7 (green 
dotted line in Figure 5). On the other hand, as long as 
QSE is not attained, the models with larger i?i yield fi¬ 
nal compositions with heavier mean mass number except 
for the model B that initially does not have heavy nuclei 
(Figure 2). 

Table 2 shows the 10 most abundant nuclei and their 
mass fractions for the model A, B, C, and D with (i?,, 
tTi)=(10, 10). In the model A, C, and D, ®®Si is most 
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Mean mass number (b) 



1 23456789 10 

Ri,7 

Mean mass number (d) 



1 23456789 10 

Ri,7 


Fig. 5. Mean mass number of the final nuclear composition in the outflow as functions of R[ and ai for (a) model A, (b) model 
B, (c) model C, and (d) model D. The model parameters are summarized in Table 1. The white solid lines and green dotted lines 
represent the criterion for the QSE establishment and the contour of Ti g = 4.7, respectively. 


Table 2. The 10 most abundant nuclei and their mass fractions for the models A, B, C, and D with {R[, cri)=(10, 10). 



Model A 

Model B 

Model C 

Model D 

1 


••• 3.7182E-01 

^He 

•• 5.1509E-01 


••• 3.9266E-01 


•• 3.8889E-01 

2 

160 

• • 2.9626E-01 

58Ni 

• • 9.0864E-02 

160 

••• 2.0150E-01 

32 s 

••• 1.7979E-01 

3 

32s 

••• 2.0439E-01 

32s 

••• 7.2190E-02 

32s 

••• 1.7598E-01 

160 

•• 1.6323E-01 

4 

36 Ar 

••• 5.3067E-02 

28Si 

• • 5.8950E-02 

54Fe 

••• 6.9657E-02 

54Fe 

•• 8.7961E-02 

5 

40Ca 

••• 4.1669E-02 

54Fe 

•• 5.7403E-02 

36 Ar 

••• 3.9364E-02 

36 Ar 

•• 4.0427E-02 

6 


••• 1.0189E-02 

55 Co 

••• 3.3212E-02 

40 Ca 

••• 2.2114E-02 

40 Ca 

••• 2.4969E-02 

7 

24Mg 

••• 7.1779E-03 

36 Ar 

••• 2.4232E-02 

55Co 

••• 1.6240E-02 

55Co 

••• 1.9738E-02 

8 

56Ni 

••• 2.8607E-03 

stni 

••• 1.9790E-02 

24Mg 

••• 1.0641E-02 

56Ni 

••• 1.4527E-02 

9 

55 Co 

••• 2.1998E-03 

160 

•• 1.8796E-02 

58Ni 

••• 8.9718E-03 

58Ni 

••• 1.3446E-02 

10 

58Ni 

••• 1.2994E-03 

56Ni 

••• 1.4453E-02 

56Ni 

••• 8.8717E-03 

4He 

••• 1.0203E-02 
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Cooling radius Rcooi [cm] 



Fig. 6. The cooling radius i?cooli above which QSE is no 

longer established, as functions of ij; and cr; for the model A. 

abundant and its mass fraction reaches ~ 40%. The mass 
fractions of and are also high in the model A, 
C, and D. In particular, the mass fraction of reaches 
~ 30% in the models A. On the other hand, in the model 
B, the most abundant nucleus is ^He the mass fraction 
of which is about ^ 52% and the mass fractions of each 
metal nucleus are less than 10%. 

So far, we fix the isotropic jet luminosity as Tj’*52=4 
and 10 referring to canonical GRB isotropic luminos¬ 
ity. However, there are GRBs belonging to less energetic 
class called low-luminosity GRBs (LLGRBs) with L!™ ^ 
1049 ej-g g-i (e.g., Liang et al. 2007). From Equation (4) 
and the equation of state, the isotropic luminosity of 
a GRB jet is related to the initial temperature through 
Cint with Equation (1) as ~ 167ri??eintr?f;i(l-|-cri)/3. 
Applying the condition for Ti to avoid the QSE estab¬ 
lishment (i.e., Ti^g < 4.7), the condition in terms of the 
isotropic luminosity is obtained as follows, 

Lf°<3.9x 10™Rf.7(l-bcri) ergs-\ (10) 

According to this condition, the metal nuclei can be in¬ 
volved in the jets of the LLGRBs even within the frame¬ 
work of standard fireball model (i.e., (Ti = 0). 

Figure 6 shows the cooling radius i?cooi above which 
the condition to establish QSE is no longer satisfied; 
.Rcooi ~ 10^“® cm depending on Ri and (Ti. Since QSE 
is not attained at r > Rcooi, the metal nuclei, that have 
not been dissociated at r ^ Rcooi, can survive after the ac¬ 
celeration. Furthermore, if the metal nuclei mix from the 
stellar mantle at r > Rcooi, as discussed in Horiuchi et al. 
(2012), they are involved in the outflow without serious 
dissociation. 

4. CONCLUSION 

In this paper, we investigated nuclear composition of 
GRB jets assuming that the jets initially possess metal nu¬ 
clei. We calculated fallback in a relativistic jet-induced ex¬ 
plosion with a two-dimensional relativistic hydrodynamics 
calculation and derived the initial composition of the jet. 
Then, we calculated the acceleration of magnetized GRB 


jets and detail nuclear reactions in the jets with the initial 
compositions. 

We found that the composition of the falling matter 
after the jet injection depends on the radius to which the 
matter falls off and thus the initial composition of the 
jet depends on the size of the central engine of the jet. 
If the size of the central engine is larger than Ri j ^ 5, 
the matter contains metal nuclei abundantly except for 
the model with Z = 0 and = 10 (model B). Model 
B involves only a small fraction of metal nuclei with the 
mean mass number of {A) < 6 even for Ri 7 ^ 5 because 
the jets consist of mass elements with high temperature 
in the presupernova star. 

We conclude that the metal nuclei can survive in the 
jet if QSE is not established. The metal nuclei are disso¬ 
ciated mostly to ^He once QSE is established. This is due 
to the high entropy of the jet being accelerated to F = 100, 
in which (a-rich freezeout takes place. Therefore, the fi¬ 
nal nuclear composition of the jet is dominated by “^He 
for the QSE-established models. The criterion for QSE 
establishment is well fitted by the contour of Ti^g = 4.7. 
The criterion leads to the condition of the isotropic jet lu¬ 
minosity Lf° as Lf° < 3.9 x 10®°R?7(1 -|-(Ti) erg s~^ and 
this is consistent with the results of Horiuchi et al. (2012) 
since the beaming correction for the jet reduces by 
the orders of ^ 2. 

The most popular model for acceleration of GRB jet is 
the fireball model, in which thermal pressure accelerates 
the jet from sub-relativistic to ultra-relativistic (Goodman 
1986; Paczyhski 1986; Meszaros 2006). In such the stan¬ 
dard fireball model, Rij > 5 is required for the metal nuclei 
to survive in the jet with ^1- If the jet has a lu¬ 
minosity of Rj ®52 ^ 4, the size of the central engine must 
be larger than 10® cm for the survival of metal nuclei. 
However, this violates a constraint on the size of central 
engine from time variability of the flux, i.e. Ri ^7 ^ 10. 
Therefore, at least, metal nuclei initially contained in the 
GRB jet with Rj ™2 ^ 4 should be destroyed and the fi¬ 
nal nuclear composition is dominated by light nuclei and 
free nucleons within the framework of the standard fireball 
model. 

On the other hand, the magnetized jet has been pro¬ 
posed to explain spectra of some GRBs. For example, 
Zhang & Pe’er (2009) suggests that GRB 080916G in¬ 
volves a magnetized jet because a thermal component, 
that should appear in the spectrum if the jet energy is 
initially dominated by thermal energy, was not detected. 
Also, Guiriec et al. (2011) reports that GRB 100724B ex¬ 
hibits typical non-thermal spectrum, called Band spec¬ 
trum, with a significant thermal component and suggests 
that a highly magnetized jet can explain the feature, which 
is quite challenging for the standard fireball model. In 
such a magnetized jet, the energy is initially possessed 
by the magnetic field, and thus the condition for QSE 
is avoidable in the jet even with high luminosities like 
Rj ®52 ^ 10 with satisfying Ri _7 < 10, although Ri .7 ^ 5 is 
required for metal nuclei to initially exist in the jet in our 
model. 





5. DISCUSSION 

In this paper, we constrain the size of the central engine 
with the time variability of the gamma-ray emission since 
it could reflect the time variability of the central engine, 
for example, in the internal shock model. However, the 
time variability of the gamma-ray emission may come from 
other factors. For examples, the order of ^ 1 s variability 
can arise from the interaction between the jet and the pro¬ 
genitor star (Morsony, Lazzati, & Begehnan 2010), or the 
time variability may be related to the emission mechanism 
itself, e.g. turbulent motion at the gamma-ray emitting 
region (e.g., Zhang & Yang 2011). In these cases, the 
size of the central engine is not necessarily constrained, 
at least, by the time variability and the larger i?i _7 > 10 
is possible. If the size of the central engine is as large as 
i?i ,7 = 50, the criterion for the survival of the metal nuclei 
is Lj®° ^ 10®^(1 -l-CTi) erg s“^, indicating that a large frac¬ 
tion of GRBs can possess the metal nuclei in the jet even 
with (Ti = 0. 

We treated the GRB jet as a steady, radial, magne¬ 
tized outflow for simplicity. However, in reality, the jet 
should propagates through the progenitor and the interac¬ 
tion between them is expected. Such an interaction could 
affect the dynamics of the jet. For examples, it is sug¬ 
gested that the initial confinement by the collapsing stel¬ 
lar material has an important role on the collimation of 
magnetically dominated jets from magnetars (Uzdensky 
& Macfadyen 2007). The cocoon, which is a shocked hot 
gas surrounding the jet, also have a role to collimate the 
jet (e.g., Morsony, Lazzati, & Begelman 2007; Mizuta & 
Aloy 2009; Bromberg et al. 2011). If the jet is collimated, 
the evolution of jet cross section E(r) may be expressed 
as S(r) oc with ^ < 2. In this case, the temperature, 
and also the density, will decrease slowly than in the ra¬ 
dial flow and the criterion, Equation (10), possibly tighten 
because the fluid will tend to keep the high temperature 
at the initial stage. 

Metzger, Giannios, & Horiuchi (2011) investigated the 
nucleosynthesis from free nucleons in the magnetically 
dominated jet in the context of the protomagnetar model 
for the central engine of GRBs (Metzger et al. 2011). The 
main difference between their work and this paper is the 
parameter range of the entropy per baryon. They esti¬ 
mated the entropy per baryon in the jet with the ana¬ 
lytic expression derived in Qian & Woosley (1996) for the 
neutrino-driven winds. They suggested that, since the 
entropy per baryon in the neutrino driven wind is suffi¬ 
ciently low, the heavy nuclei beyond the iron peak, i.e., 
A ^ 56 where A is a mass nubmer, can be synthesized in 
the jet. On the other hand, we do not consider such suf¬ 
ficiently low entropy environment for synthesis of metal 
nuclei from free nucleons and we focus only on the sur¬ 
vival of the metal nuclei in this paper. Future studies on 
the central engine of GRBs will reveal whether such low 
entropy enviornment is realized or not. 

It is known that GRBs are one of the candidates for the 
origin of ultra high energy cosmic rays (UHEGRs) (e.g., 
Waxman 1995; Vietri 1995) and LLGRBs are also possible 


candidates for the origin of UHEGRs (e.g., Murase et al. 
2006; Gupta & Zhang 2007). The acceleration of metal 
nuclei up to ultra-high energies could be possible in both 
usual high-luminosity GRBs and LLGRBs (Wang et al. 
2008; Murase et al. 2008). Abraham et al. (2010) reported 
that the observed UHEGRs are dominated by heavy nuclei 
at high energies, i.e., E > 10^® eV (see, however, Abbassi 
et al. 2010). Hooper & Taylor (2010) suggested that the 
result of Abraham et al. (2010) can be quantitatively re¬ 
produced only if UHECR composition at the accelerating 
site mainly consist of intermediate mass nuclei such as ni¬ 
trogen together with a considerable fraction of heavy nu¬ 
clei such as iron. Therefore, GRBs with magnetized jets 
or LLGRBs have a potential to be the origin of UHEGRs 
in terms of, at least, the nuclear composition because in¬ 
termediate mass nuclei, e.g. ^®Si, are very abundant in 
the jet. However, in order to explore whether the model 
considered in this paper can quantitatively explain the 
UHEGR observations, acceleration process, escape from 
the source, and the propagation of the UHE nuclei from 
the source to the Earth also have to be considered. Such 
calculations are beyond the scope of this paper and will 
be studied elsewhere. 
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